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Abstract 

The study of quantum Coulomb systems at equilibrium is important for understanding 
properties of matter in many physical situations. Screening, recombination and van der 
Waals forces are basic phenomena which result from the interplay of Coulomb interactions, 
collective effects and quantum mechanics. Those phenomena are introduced in the first 
part of this lecture, through various physical examples. Their treatment within mean- 
field theories and phenomenological approaches is also exposed, while related predictions 
are discussed. This sheds light on fundamental issues, which must be analyzed without 
any a priori approximations or modelizations. The second part of this lecture is pre- 
cisely devoted to the presentation of various exact results for the quantum proton-electron 
hydrogen plasma. Such results are derived within the Screened Cluster Representation, 
which is constructed by combining the path integral representation of the Coulomb gas 
with Mayer-like diagrammatical techniques. They illustrate the breakdown of Debye ex- 
ponential screening by quantum fluctuations, as well as the emergence of familiar chemical 
species in suitable low-temperature and low-density limits. Also, the amplitude of van der 
Waals forces is shown to be reduced by free charges. 

PACS numbers: 05.30.-d, 05.70.Ce, 52.25.Kn 
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I. INTRODUCTION 



Under standard Earth conditions, and also in many astrophysical situations, the 
properties of matter result from the interplay between non-relativistic quantum me- 
chanics and Coulomb interactions. All relativistic effects, as well as the other non- 
electromagnetic interactions can be safely omitted. A description of matter in terms 
of quantum nuclei and quantum electrons interacting via the Coulomb potential is 
then sufficient. In that context, the derivation of exact results for equilibrium prop- 
erties of quantum Coulomb systems is of crucial importance. 

The study of quantum Coulomb systems within statistical mechanics, requires to 
face several difficult problems related to both short- and long-range specificities of 
Coulomb potential, as well as to quantum mechanics itself, namely screening and 
recombination. In the first introductory part of the present lecture, those problems 
are successively adressed for various physical examples where free charges interacting 
via the Coulomb potential are present. First, it is instructive to consider screen- 
ing in classical systems where quantum effects can be omitted. According to the 
pioneering Debye-like mean-field theories, Coulomb interactions are exponentially 
screened at the classical level, a prediction confirmed by a large variety of rigorous 
proofs and exact results. In the quantum case, similar mean-field theories also pre- 
dict an exponential decay of equilibrium particle correlations. However, they do not 
account for a very fundamental feature of quantum mechanics, namely the intrusion 
of dynamical effects at equilibrium. Hence, there exist various arguments which sug- 
gest a breakdown of exponential screening in quantum systems. In the literature, 
recombination has been mainly dealt with in the framework of the chemical picture. 
The corresponding phenomenological approaches are based on ad hoc modelizations 
for preformed entities and their interactions. 

As illustrated by the various considerations exposed in the introductory part of 
this lecture, screening, recombination and van der Waals forces result from entangled 
mechanisms combining Coulomb interactions and quantum mechanics. Although, 
phenomenological treatments of such phenomena have been widely developped, there 
still remain subtle effects, a deeper understanding of which requires a more funda- 
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mental analysis in the framework of quantum Coulomb systems, namely the so-called 
physical picture. The aim of the second part of this lecture is to present several exact 
results relative to fundamental issues about screening and recombination. For the 
sake of simplicity and pedagogy we consider the quantum hydrogen plasma made 
with point protons and point electrons interacting via the Coulomb potential. First, 
we present rigorous proofs about thermodynamical stability on the one hand, and 
the atomic limit on the other hand. Then, we introduce the Feynman-Kac path 
integral representation, which turns out to be a quite efficient tool for our purpose. 
Its application to a many-body quantum system provides an equivalent classical sys- 
tem made with extended objects called loops. Standard Mayer-like diagrammatical 
series for the gas of loops, are exactly transformed into the Screened Cluster Repre- 
sentation for equilibrium quantities of the quantum system. That transformation, 
based on suitable resummations and reorganizations, accounts simultaneously for 
both screening and recombination. Quantum fluctuations, which play a crucial role 
in screening, are merely embedded in loop shapes. Recombination, which cannot 
be treated perturbatively in the Coulomb potential, is automatically ensured by 
the presence of Boltzmann factors associated with loop interactions. The Screened 
Cluster Representation is applied to hydrogen in the Saha regime, defined by both 
low temperatures and low densities, where it behaves as a partially ionized atomic 
gas. An exact asymptotic expansion for the equation of state is constructed, be- 
yond familiar Saha theory. It sheds light on a suitable first-principles account of 
contributions of usual chemical species, without any a priori modelizations. Also, 
particle correlations are found to decay algebraically. Thus, ionized protons and 
ionized electrons only reduce the amplitude of van der Waals interactions between 
atoms. As a conclusion, we summarize the main answers to the above issues inspired 
by those exact results. 
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II. COULOMB INTERACTIONS AND QUANTUM MECHANICS AT 
WORK 



A. Examples and specific features 

A description of matter in terms of Coulomb systems, also called plasmas, is 
necessary as soon as free charges are present. There is a large variety of physical 
examples where such situations occur. For instance, electrolytes involve ionic species 
obtained by dissolution of salts into water. Also, in metals, the electrons of the 
conduction band freely move across the samples. Other examples can be found in 
astrophysical situations, where high temperatures or high pressures favor ionization 
of matter, in general inside the cores of compact objects or stars. 

Let us consider a plasma made with point charges. Two charges g; and qj located 
at Ti and rj interact via the instantaneous Coulomb potential 

«c(ri,r i ) = r ^- r , (HI) 

while the full interaction potential of the system reduces to the sum of pairwise inter- 
actions ( IH.ip . Here, we dot not taken into account retardation effects or magnetic 



forces, and we discard the coupling of charges to elecromagnetic radiation. That 
purely Coulombic description is sufficient as far as the average speed of charges is 
small compared to the speed of light c. Roughly speaking, for a classical charge with 
mass m, this requires that the thermal energy fc^T is small compared to the rest 
energy mc 2 This implies T <C lO 10 -^ for electrons and T <C 10 13 i^ for protons. 
Such conditions are fullfilled in many physical systems! 

Contrarily to interactions between neutral entities, Coulomb potential is long 
ranged, namely it is not integrable at large distances 

dr- = oo, (II.2) 

r>D r 

where D is some irrelevant cut-off length. That divergence might pollute thermo- 
dynamical quantities, and usual extensivity properties might be lost. The underly- 
ing physical picture is that charges with the same sign tend to repell together too 
strongly at large distances, so the system might explode. 



In addition to above long-range problems, the 1/r-nature of the Coulomb poten- 
tial causes also some trouble at short distances. If we consider two classical opposite 
point charges ±q separated by a distance r, the corresponding Boltzmann factor is 
not integrable at r = 0, i.e. 



with (3 = 1/ (ksT). Thus, a classical system with positive and negative point charges 
collapses. 

According to above considerations, the study of Coulomb systems in the frame- 
work of statistical mechanics has to face two central difficulties related to the be- 
haviours of the Coulomb potential at respectively large (r — > oo) and short (r — > 0) 
distances. In the following, we will describe the fundamental mechanisms which cure 
the corresponding singularities, namely screening and recombination. Screening is 
a collective effect which prevents explosion at large distances. Recombination is a 
purely quantum mechanical effect, where the uncertainty principle smears out the 
divergency of the Coulomb potential at r = 0, so collapse is avoided. Those mecha- 
nisms are adressed successively, by considering physical systems for which simplified 
models and/or phenomenological approaches can be introduced. 

First, in Section III BL we present Debye theory for a classical electrolyte. That 
mean-field approach predicts an exponential screening of Coulomb interactions, 
which is confirmed by rigorous results. A similar mean-field theory has been applied 
to the quantum electron gas, as described in Section III CI Again, an exponential 
screening is found, but that prediction is shown to be quite doubtful according to 
the unavoidable intrusion of dynamical effects in equilibrium quantities of quantum 
systems. Recombination and van der Waals forces are considered in Section III Dl 
through the example of hydrogen in the Sun, the state of which changes from a fully 
ionized plasma in the core to an atomic gas in the photosphere. We present simple 
arguments, as well as some standard phenomenological considerations, which intro- 
duce the more sophisticated analysis presented in the second part of this lecture. 
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B. Screening in classical systems 



1. Primitive model for an electrolyte 

Let us consider an ordinary solution of sodium chloride. For a concentration 
Co = 0.1 moles/1 at room temperature T = 300 K, the salt is entirely dissociated, 
so the solution reduces to a mixture of ions Na + , Cl~ and water molecules H 2 
(ions H 3 + and OH~ are omitted since their concentration is negligible compared 
to C )- The mean interionic distance, a = (3/ (An p)) 1 ^ with the common ionic 
number densities p = p + = p_, is large compared to both the typical size d of ions 
and the mean distance between water molecules. Therefore, a large number of water 
molecules surrounds each ion, and water may be reasonably replaced by a continuous 
medium with dielectric constant e w . Also the de Broglie thermal wavelengths A± = 
(/3h 2 /M±) 1 / 2 are small compared to d, so ions can be treated classically. 

According to the above hierarchy of length scales, the present electrolytic solu- 
tion is well described by a classical two-component plasma made with hard spheres 
carrying charges q± = ±e, also called the restrictive primitive model in the litera- 
ture. Two ions of species a — ± and 7 = ± separated by a distance r interact via 
the potential 



The hard-core part of the potential is a phenomenological description of the short- 
range repulsion between the electronic clouds of the two considered ions. Of course, 
within that modelization, no collapse between oppositely charged ions occur. Also, 
notice that the bare Coulomb potential is renormalized via the standard factor l/e w 
which accounts for the underlying presence of a continuous dielectric medium. 

In order to adress the question of screening, we introduce the equilibrium charge 
density C a (r) surrounding an ion with species a fixed at the origin, which reads 



t( a7 (r) = +00 , r < d a 
u ay (r) = , r > d, 



(II.4) 




(II.5) 
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FIG. 1: A typical ionic configuration C where two cations Na + (0 = 7 = +) are fixed at 
and r respectively. 

(21 

where p a -y is the two-point equilibrium distribution of species a and 7. Screening 
can be easily understood within a mean-field calculation of C a (r). The basic ideas 
sustaining nrean-fie.d approach have been first .ntrodnced by Gony Q and Chap- 
man |2[ fot the study of electrical double layers near charged electrodes. They have 
been extended to the calculation of polarization clouds in the bulk phase by Debye 
Here, we present the corresponding arguments which provide the 



and Huckel 



a. 



mean-field form of C a (r). 



2. Debye theory 

Let C be a given spatial configuration of ions where one ion a is fixed at the 
origin and one ion 7 is fixed at r. In Fig. HJ we draw a typical configuration which 
illustrates that ion a attracts oppositely charged ions and repells the other ones. Ion 
7 feels the electrostatic potential <p(r\C) created at r by all the remaining ions, in 
particular those which are far apart, because the Coulomb potential is long ranged. 
Since f{r\C) is a sum of a large number of terms, fluctuations of the positions of 
non-fixed ions should slightly modify its value. That statement is inspired by the 
law of large numbers, and it can been applied to other many-body problems where 
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long range interactions are present [751 ] . Thus, it is tempting to replace <p{v\C) by its 
average value <p a (r), which is nothing but the electrostatic potential created by ion 
a plus its polarization cloud with charge distribution C a (r). Within that mean-field 
treatment, the density of ions 7 at r is merely given by the Boltzmann distribution 
of an ideal gas in the external potential g 7 <p a (r), i.e. 

pS(0,r) 



Pa 



p 7 exp(-/3g 7 y3 a (r)) . (II. 6) 



At this level, we have completely omitted short-range effects arising from the 
hard-core part of ion-ion potential (111.4j) . Under diluted conditions for which d <C a, 
this is legitimate since for most configurations C, ion 7 only feels the Coulomb parts 
of its interactions with other ions. At the same time, position fluctuations of neigh- 
bours of fixed ion 7 can be neglected, only if they generate variations of (p(r\C) which 
do not exceed the thermal energy ksT. This implies that the corresponding typi- 
cal value e 2 /(e w a) of Coulomb interactions is small compared to ksT, a condition 
indeed fulfilled at low densities. According to previous considerations, the lineariza- 
tion of the Boltzmann factor ( 1II.6j) is consistent with the mean-field treatment, and 



it provides 

C a (r) = -/3^g> 7 ^(r), (II.7) 

7 

where we have used the neutrality condition 

7 

The insertion of the linearized expression ( 111.71) of C a (r) into Poisson equation for 
(f a (r), shows that <p a (r) = q a (j)-£,(r), where </>d is the solution of 

4-7T 

(-A + 4) Mr) = - 6(t) (U.9) 

6w 



with 



and boundary conditions 



4 = -/3E*7 = -/3e 2 p (11.10) 

7 



6 D (r) ->■ when r — >• 00 . (H.ll' 
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Such boundary conditions amount to impose that the ionic densities induced at r 
tend to the bulk homogeneous value p when r — > oo, or in other words to make the 
quite plausible assumption that polarization effects decay at large distances. 



The Helmholtz equation (1II.9j) together with boundary conditions is easily solved 



in Fourier space, and 0D( r ) is then readily obtained by applying Cauchy's theorem 



see e.g. Ref. 6] for a detailed presentation of that standard calculation), 

M r) = Mr) = • (H.12) 

Thus, the electrostatic potential q a 4>n{ r ) inside the electrolytic solution decays ex- 
ponentially faster than the potential q a / (e w r) created by a single ion a immersed in 
water. In the electrolyte, such an ion is screened by its induced polarization cloud 
over the so-called Debye length Ad = k^ 1 . 

3. Reliability and limits of mean-field predictions 

It is instructive to determine what are the physical conditions under which the 
mean-field treatment is reliable at a quantitative level. The main argument for 
neglecting fluctuations relies on the large value of the number of ions which con- 
tribute to the electrostatic potential at r. In average, those ions are contained in a 
sphere with radius r as a consequence of Gauss theorem, so their number is of order 
8irpr 3 /3. That number must be large for r of order the typical length scale Ad- This 
implies okd <C 1, or equivalently the weak-coupling condition e 2 / (e^aksT) <C 1. 
Within that condition, the linearization of Boltzmann factor flll.6j) is also justified. 



If that condition is fulfilled at high temperatures or low densities, hard-core effects 
can be neglected only if d <C a. The previous mean-field approach is then expected 
to be valid, at least at a quantitative level, at sufficiently high temperatures and 
low densities. Notice that in the present example, the weak-coupling condition is 
ensured thanks to the large value of the water dielectric constant, e w ~ 80, which 
drastically reduces the strength of electrostatic interactions. As a consequence, De- 
bye theory works reasonably well for most electrolytic solutions of ordinary salts 
under standard conditions. 
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Beyond its practical interest for a large variety of physical systems, Debye the- 
ory suggests a fundamental result, namely the exponential decay of correlations in 
classical charged systems, at least in fluid phases. In fact, that remarkable property 
has been proved for various systems {7 ^] at sufficiently high temperatures and suffi- 
ciently low densities, i.e. in thermodynamical regimes where mean- field approach is 
expected to work. The very difficult part of those proofs relies on the complete con- 
trol of the contributions of all the effects omitted in Debye theory. It turns out that 
such effects do not destroy the exponential clustering predicted by mean-field. In 
other words, the essence of classical exponential screening is captured by mean-field 
approach, which enlights its collective nature. Also, the harmonicity of the Coulomb 
potential is a key ingredient, as illustrated by an analysis of the equilibrium BGY 



hierarchy lOj. In general, the fast decay of particle correlations is proved to be 
related to a perfect arrangment of polarization clouds, as exemplified through mul- 
tipole sum rules jll]. Here, the total charge of the polarization cloud surrounding 
a must exactly cancel out its charge, namely 

drC a (r) = -q a . (11.13) 

That monopole sum rule is indeed satisfied by the mean-field expression 

2 exp(-/c D r) 

-q a K D — (11.14) 

of C a (r). Of course, the existence of positive and negative charges is crucial for 
screening mean-field expression (III. 14ft diverges 

when r — » 0. Such divergency is unphysical, and it proceeds from the linearization 
of Boltzmann factor ( 111.61) which is never valid at short distances. However, that 
spurious short-range singularity does not affect integrated quantities, like the total 
charge ( III. 13ft or thermodynamical quantities. 



At a quantitative level, corrections to Debye theory can be derived within Abe- 
Meeron diagrammatical expansions [ijj]. Debye theory is merely recovered by 



keeping only the first graph in those series. Each of the remaining graphs provides an 
exponentially decaying contribution to particle correlations, in agreement with above 
proofs. The resulting screening length As(p, T) differs from its Debye expression Ad, 
and it might be approximately estimated by selecting suitable classes of Abe-Meeron 
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graphs. Notice that short-range repulsion contributes to Xs(p,T) in a rather subtle 
way, so a reliable description of the variations of Xs(p,T) when temperature is 
decreased and/or density is increased is a challenging problem. At sufficiently low 
temperatures and/or high densities, the occurence of phase transitions might lead 
to divergencies of \$(p, T). For instance, exponential screening might be lost in a 
crystalline phase, or at the critical point of the liquid-gas transition as suggested in 



Ref. 



16|- 



Eventually, let us conclude that section by a very important remark about dy- 
namical aspects. In a real sample, ions continuously move, and because of their finite 
inertia, instantaneous perfect arrangments of the polarization clouds cannot be en- 
sured. However, according to the somewhat magical recipe of statistical mechanics, 
at equilibrium the time-average of any quantity can be equivalently computed as 
a phase-space average with the Gibbs measure. Thanks to the remarkable fac- 
torization of that measure into the Maxwellian distribution of momenta times the 
Boltzmann factor associated with the interaction potential, dynamical features are 
washed out by momenta integration. Then, the calculation of equilibrium correla- 
tions in position space becomes a purely static problem, where only configurations 
of particle positions must be considered without caring about particle velocities. In 
other words, dynamical fluctuations no longer intervene, so above considerations 
about particle inertia are not relevant. However, in a quantum system where above 
factorization of the Gibbs measure is no longer valid, we can anticipate that dynam- 
ical effects should contribute to equilibrium correlations. 



C. Screening in quantum systems 

1. Jellium model for conduction electrons 

In ordinary metals, like Copper for instance, there exists a finite density p of 
ionized electrons which ensure electrical conduction. Those electrons freely move 
among the remaining ions, which can be considered as fixed at their lattice sites. A 
further simplification is to replace the periodic ionic charge density by the constant 
ep. This provides a model for conduction electrons, called either jellium or one- 
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component plasma, where such electrons are immersed in an uniform and rigid 
neutralizing background. 

Under standard conditions and for simple metals with one or two ionized electrons 
per atom, the typical values for temperature and density are T = 300 K and p = 
3.10 28 m~ 3 . The thermal de Broglie wavelength of electrons A e = ((3h 2 /m c ) is then 
large compared to the mean interelectronic distance a = (3/(47rp)) 1//3 , so electrons 
must be described by quantum mechanics. Equivalently, thermal energy /c#T is small 
compared to Fermi energy Ef = h 2 k F /(2m e ) with Fermi wavenumber kp = (37r 2 p) 1 / 3 , 
so electrons are strongly degenerate and their typical kinetic energy is of order £p. 
Since e-p is small compared to the rest energy m e c 2 , relativistic effects can be omitted. 
Therefore, we have to consider the non-relativistic quantum version of jellium, where 
two electrons separated by a distance r interact via the instantaneous Coulomb 
potential 

Mee(r) = C - • (11.15) 

In addition, each electron is submitted to the electrostatic potential created by the 
ionic background, while the constant background self-electrostatic energy is also 
taken into account in the full interaction potential of the system. Electrons obey 
to Fermi statistics since they carry an half-integer spin. Notice that the Coulomb 
Hamiltonian of the present model do not depend on electron spins. 



2. Thomas- Fermi theory 



Similarly to the classical case described above, screening properties can be merely 
illustrated through a mean-field calculation of the equilibrium charge density C e (r) 
surrounding an electron fixed at the origin, 



Ce(r) 



, pff(0,r) 
-< I p 



(11.16) 



In definition flII.16j) . p^ is the two-point equilibrium distribution of electrons, while 
the substracted term accounts for the ionic background density. The so-called 
Thomas- Fermi theory 17|, ll8j was introduced for describing the electronic structure 
of atoms, and it can be viewed as the source of further density functional methods. 
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Its application to the study of screening properties [19| can be rephrased similarly to 
the mean- field approach for the classical case. Namely, we assume that the electrons 
close to a given point r constitute an ideal gas submitted to the external one-body 
potential — e<p e (r), where <p e ( r ) is nothing but the electrostatic potential created by 
the electron fixed at the origin plus its polarization cloud with charge distribution 
C e (r). Now, according to the quantum nature of that ideal gas, its density is related 
to the external potential via the Fermi-Dirac distribution instead of the Boltzmann 
law, i.e. 

A0,r) = J2_ f 1 

p ~ (2vr)3 J exp|)0(e(k) - e^ c (r) -//)] + 1 ' ' * ' 

In Fermi-Dirac expression (III.17j) . e(k) = h 2 k 2 /(2m e ) is the non-relativistic kinetic 
energy of a plane wave with wavenumber k, while the chemical potential p, taken 
constant across the system, is naturally determined by the bulk condition 

"=(2^/ dk exp[/3( E (k!- A )] + l- (IU8) 

Assuming that polarization effects vanish at large distances, we set pee^(0,r) — > p 2 
when r — > oo. The corresponding limit form of equation (III. 171) then implies the 
boundary condition (p e (r) — > 0. 

As argued for the Boltzmann factor involved in Debye theory, the linearization of 
Fermi-Dirac distribution (III. 171) with respect to <p e (r), is consistent with the present 
mean- field treatment, which also requires that Coulomb interactions are small per- 
turbations. If we set v ? e( r ) = ~ e0TF( r ) ; we then obtain an Hemholtz equation for 
0TF( r ) which is identical to its classical counterpart f ]II.9j) for 0o( r ), except for the 



replacements e w — > 1 and — > htf where the Thomas-Fermi wavenumber ktf is 
defined by 



2 d P. 



4 F =4vre 2 ^(/3,/i). (11.19) 
The resulting mean- field expression for charge distribution C c (r) reads 

_ e4 exp(-. TF r) 
47rr 

Thus, the mean-field approach again predicts an exponential screening. Now the 
screening length Atf takes into account degeneracy effects controlled by dimension- 
less parameter ep/ksT. In the limit of weak degeneracy, which can be obtained by 
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setting 11 — > — oo at fixed (3, the density vanishes as p ~ 2exp(/3/z)/(27rA e ) 3//2 , so 
Atf does reduce to its classical Debye form (4ii(3e 2 p)~ 1 / 2 . In the opposite strong- 
degeneracy limit, obtained by setting p — > +00 at fixed (3, the density diverges as 
p ~ (2m c /i) 3/2 /(37r 2 ) and we find 

5. Expected validity regime 

As quoted above, the validity of Thomas-Fermi theory requires weak-coupling 
conditions. Of particular interest is the regime of strong degeneracy, where the 
relevant kinetic energy is the Fermi energy e-p. The corresponding weak-coupling 
condition reads e 2 /{ae-p) <C 1. That strongly-degenerate weakly-coupled regime is 
reached at sufficiently high densities for a given temperature. Indeed, when p — > +00 
at fixed 0, Fermi energy obviously becomes larger than classical thermal energy, 
while Coulomb energy of order p 1//3 grows slower than kinetic Fermi energy of order 
p 2 / 3 . In that high-density limit, notice that Atf, given by expression ( III. 211) . becomes 
large compared to a, so the number of electrons inside the screening sphere with 
radius Atf is indeed large. 

In order to derive corrections to the mean-field approach, it is necessary to use the 



formalism of many-body perturbation theory 20]. In that framework, the mean-field 
ideas are recasted into the celebrated Random Phase Approximation (RPA), which is 
close to Thomas-Fermi theory. In particular, RPA also predicts an exponential decay 
of charge correlations. Systematic corrections to RPA can be derived by taking into 
account suitable Feynman graphs. This leads to asymptotic high-density expansions 
of thermodynamical quantities, like the pressure or the internal energy, beyond their 
ideal strongly-degenerate forms. 

For the strongly-degenerate electron gas described above, the Coulomb energy 
is not small compared to the Fermi energy. Therefore the system is not weakly 
coupled, and the application of Thomas-Fermi theory is rather questionable. In 
particular, we find that Atf — 0.1 nm is smaller than a ~ 0.2 nm, in violation of 



the mean-field validity condition Atf ^> a \77\ 
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4- Doubts about mean-field predictions 



If mean-field approaches like RPA may provide the leading behaviour of inte- 
grated quantities involving C e (r) in the high-density limit p — » +00 at fixed (3, 
their reliability for C e (r) itself is not uniform with respect to distance r. At short 
distances, the divergency of Thomas- Fermi expression ( 1H.20|) for C c (r) is obviously 
unphysical, like for its classical counterpart in the framework of Debye theory. This 
is due first to the strong repulsion embedded in </3 e (r) when r — > which cannot be 
treated perturbatively. Also TF theory does not account for the effective repulsion 
induced by Pauli principle when both electrons located at and r have the same 
spin orientations. At large distances, and contrarily to the classical case, mean-field 
predictions are also doubtful, as argued below. 



In writing the Fermi-Dirac distribution f ]II.17j) . a very important implicit assump 



tion is made. Indeed, potential y2 e (r) is treated as a constant, so the energy levels 
are merely shifted by that constant from their corresponding purely kinetic expres- 
sion e(k) = h 2 k 2 /(2m e ). This amounts to neglect the spatial variations of </9 e (r), 
or in other words to omit diffraction effects arising from the non-commutativity 
of the kinetic and potential parts of Hamiltonian —H 2 A/(2m e ) + <p e (r). Then, we 
are left with a purely configurational static problem, like in the classical case, and 
not surprisingly this leads also to an exponential screening. However, that result 
is rather questionable since an intrinsic specificity of quantum mechanics has been 
erased from the start! 

The occurence of quantum diffraction in equilibrium statistical weights is a signa- 
ture of dynamical effects, which are consequently still present at equilibrium. Now, 
and contrarily to the classical case, dynamical fluctuations combined to inertia effects 
should prevent instantaneous perfect arrangments of equilibrium screening clouds. 
Thus, the presence of partially screened multipolar interactions can be anticipated, 



and this might lead to a breakdown of exponential screening in equili 
tions. That quite plausible scenario has been first conjectured in Ref. 



Drium correla- 



22], according 



to the presence of algebraic tails in some imaginary-time Green functions. Notice 



that an analogous scenario occurs for classical time-displaced correlations 



23) 
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D. Recombination and van der Waals interactions 

1. Hydrogen in the Sun 

The sun is mainly made of hydrogen and helium, while heavier elements like oxy- 
gen, carbon, iron,... are less abundant and contribute to a small fraction of its total 
mass. Here, we will focus our study on hydrogen, and we will discard its interac- 
tions with all the other elements. Also, we do not consider complicated phenomena 
at work inside the Sun, like radiation transfer or convection for instance. For our 
purpose, it is sufficient to assume that some local thermodynamical equilibrium is 
reached at any place, and then the central question is to determine the correspond- 
ing chemical composition of an hydrogen plasma in terms of ionized and recombined 
entities. 

In the core of the Sun, temperature T is of the order of 10 6 K and density p 
is of order 1 g/cc. Since ksT ~ 90 eV is much larger than the Rydberg energy 
of order 13.6 eV, atoms, molecules, as well as other recombined entities are then 
fully ionized. The hydrogen gas is then mainly composed of free protons and of free 
electrons. In intermediate layers, temperature and density decrease as the distance 
to the core increases. This favors recombination of protons and electrons into atoms, 
and hydrogen then behaves as a partially ionized atomic gas. In the photosphere, T 
is of order 6.10 3 K while p is of order 10" 7 g/cc. Under such cool (fcsT ~ 0.5 eV much 
smaller than Rydberg energy) and diluted (the mass density of Earth atmosphere 
is of order 10~ 3 g/cc) conditions, the full proton-electron recombination is achieved, 
and hydrogen reduces to an atomic gas. That recombination process is schematically 
represented in Fig. [2j 

A proper description of above recombination process within statistical mechan- 
ics requires to account for both quantum mechanics and Coulomb interactions as 
detailed in the next Section. Here, we list several natural problems, and we present 
simple preliminary arguments which are particularly useful for a better understand- 
ing of more sophisticated analysis carried out in the framework of the many-body 
problem. 
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FIG. 2: A schematic view of the Sun. Signs + and — represent ionized protons and ionized 
electrons respectively, while circles represent hydrogen atoms. 

2. How to define an hydrogen atom? 

In quantum mechanics, all particle positions are entangled in a global wavefunc- 
tion describing the whole system. If we consider an eigenfunction of the total Hamil- 
tonian, no two-body wavefunction describing an hydrogen atom can be factored out. 
In that context, an unambiguous definition of such an atom is only possible in the 
double limit p — > and T — > 0. Indeed, the zero-density limit ensures that once a 
proton has married an electron, all other neighbours are very far apart. Also, the 
zero-temperature limit guarantees that the atom stays in its groudstate and no ther- 
mal ionization occurs. Under such infinitely diluted and cold conditions, it can be 
reasonably expected that a single hydrogen atom in the vacuum naturally emerges. 



3. How protons and electrons recombine into atoms? 



The formation of atoms in a gas of protons and electrons has been studied long 
ago by Saha [24], in the framework of the chemical picture. He considered protons, 
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electrons and atoms as three independent species with their own chemical potentials. 
If all interactions are neglected and if each species is viewed as a classical ideal gas, 
the mass action law associated with the recombination equation p + e =f± H, reads 
in terms of species ideal densities 

Pp ld) Pe d) _ exp(/3E H ) m 99 ^ 

where E-r is the groundstate energy of an hydrogen atom, En = —me 4 /(2ft 2 ), with 
reduced mass m = (m p m e )/(m p + m e ) and associated de Broglie wavelength A pe = 
((3ft 2 /m) 1 / 2 . The temperature-dependent reaction constant in the r.h.s. of mass 
action law (III. 22j) . goes to zero when T vanishes, while it explodes when T diverges. 
Therefore, and as expected, low temperatures favor atomic recombination, while 
high temperatures favor ionization. 

At given temperature and density, the chemical composition of the hydrogen gas 
is determined by adding to mass action law fill. 22 j) the neutrality rule, p p ld ^ = pe d \ 
as well as the total density constraint, p p d ^ + = p. A straightforward calculation 
then provides 



p p id) = Pf ) = P 



* (yr+v^-i) (H.23) 



and 



(id) * 
Pat = P 



(l + p/p* - y/l + 2p/p*) (11.24) 
with temperature-dependent density 

* exp(/3E H ) m v 

" = ' (II ' 25) 

According to those simple formulae, hydrogen becomes fully ionized for p < /)*, 
while it reduces to an atomic gas for p ^> p* . This explains the behaviour of 
hydrogen inside the Sun, since the corresponding ratio p/p* varies from 10~ 3 in the 
core to 10 7 in the photosphere. 



4- Why atoms and not molecules? 

In Earth atmosphere, hydrogen is a molecular gas and no atoms are present. 
This is oftenly explained by invoking that molecule H2 is more stable than atom H, 
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FIG. 3: Phase diagram of hydrogen at low densities. Regions where either ionized charges, 
atoms or molecules prevail, are separated by dashed {p*{f3)) and dash-dotted lines, in the 
vicinity of which the respective dissociation/recombination processes take place. 

namely _Eh 2 < 2 En, where E^ 2 is the molecule groundstate energy. However, we 
stress that such an argument strictly applies at zero temperature, where only energy 
U* intervene. At finite temperatures, T > 0, entropy S* must be also considered 
in order to determine the chemical composition of the system. For a given number 
N — N p — N e oi protons and electrons enclosed in a box with volume A at a given 
temperature T, the entropy of the atomic gas is larger than that of its molecular 
counterpart, because the atom number iVn = N/2 is twice the molecule number 
A^h 2 = N/4. Thus, in the minimization of the thermodynamic potential (U* — TS*) 
at fixed T, A and N, entropy contribution favors the formation of atoms [78| . That 
effect becomes more important as T increases or p = N/A decreases. In Earth 
atmosphere, temperature is rather low and density is rather high so only molecules 
are formed. In the Sun, temperature is sufficiently large and density sufficiently low 
so atoms prevail. The phase diagram of hydrogen at low densities is schematically 
drawn in Fig. El 
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5. What about the internal atomic partition function? 

At finite temperatures, not only the atomic groundstate contributes to thermo- 
dynamic quantities, but also excited states. When T vanishes, an atom can be 
considered as frozen in its groundstate because of the finite gap between and the 
first excited level Eh/4. The corresponding contribution to equilibrium quantities 
merely involves Boltzmann factor exp(— /3Eu). At finite temperatures, contributions 
of excited states must also be taken into account. For that purpose, it is natural to 
introduce the internal atomic partition function in the vacuum 

oo 

z (vac) = ^ n 2 exp{-f3E n /n 2 ) . (11.26) 

n=l 

In that definition, factor 4 arises from proton and electron spins degeneracy, factor 
n 2 accounts for orbital degeneracy of the level n with energy E^/n 2 , and summation 
extends over all boundstates up to n = oo. 

Obviously Z H vac ' ) diverges because of the contributions of Rydberg states with n — > 
oo. In the framework of the many-body problem, that divergency is usually expected 
to be cured thanks to screening effects, according to the following rough argument. 
At finite T and p, there is a finite amount of free charges, so the proton-electron 
potential inside an hydrogen atom decays faster than 1/r at distances r larger than 
screening length As- Thus, there is no accumulation of atomic boundstates near 
the continuum edge at E — 0, and all Rydberg states with n sufficiently large 
are suppressed. Since the spatial extension of such states is of order n 2 a-Q with 
Bohr radius Ob = h 2 /(me 2 ), the infinite sum ( III. 261) has to be truncated up to 
Ps ~ (As/«b) 1//2 - This provides a non-divergent atomic partition function, 

PS 

^(scr) = 4 ^ n 2 exp{-f3E u /n 2 ) , (11.27) 

n=l 

which accounts for screening effects. 

6. Are van der Waals interactions screened? 

As guessed by van der Waals in his thesis, atoms or molecules attract themselves 
at large distances. For atoms with no permanent dipole, that interaction is due to 
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quantum fluctuations as first shown by London [25( . Let us sketch the corresponding 
calculation for two hydrogen atoms in their groundstate. Each atom i (i = 1,2) is 
described by an Hamiltonian H$ , and the full Hamiltonian of both atoms reads 

H = + Kt,at + #1? (11.28) 

where V^t.at is the electrostatic interaction between the two atoms. For large sep- 
arations R between such atoms, namely R ^> Ob, that potential becomes dipolar 
because of atom charge neutrality, i.e. 

Kt.at ~ (Pi • V) (p 2 • V) - when R^oo, (11.29) 

where is the instantaneous electrical dipole carried by atom i. At such distances, 
when determining the groundstate energy of the whole two-atoms system, V at ,at 
can be treated as a small perturbation with respect to + . The first- 
order correction vanishes because the quantum average of p« in the unperturbed 
groundstate is zero. The second-order correction can be reinterpreted as an atom- 



atom effective potential which reads 
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u v MR) = ^ • (n.30) 

In that expression, CVaw is a pure numerical constant which is positive, so van 
der Waals potential U v $w(R) is attractive. That potential appears to be generated 
by quantum fluctuations of the instantaneous l/i? 3 -dipolar interaction between two 
atoms : this explains the l/i? 6 -decay of tZvdw(-R) at large distances, which merely 
follows by taking the square of 1/R 3 . Eventually, we check a posteriori that use of 
perturbation theory for deriving formula ( 1II.30j) is indeed justified if R 3> Ob, since 
\U y< iw(R)\/\Eh\ is proportional to (ob/-R) 6 - 

Above derivation is performed for two isolated atoms at zero temperature in the 
vacuum. In the framework of the many-body system at finite T and p, besides the 
above difficulties in defining properly atoms, an interesting question concerns the 
effects of free charges on van der Waals interactions. Since the instantaneous atom- 
atom potential V ati at is of purely Coulombic origin, screening mechanisms should 
modify U v dw(R) at distances R larger than screening length As. A very crude and 
naive application of the classical recipe would lead to an exponential decay of the 
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effective atom-atom interactions. Nevertheless, and as it has been anticipated in 
the construction of the phenomenological Derjaguin-Landau-Verwey-Overbeek the- 
ory for mixtures of charged and neutral fluids 27H29I]. screening of van der Waals 
interactions should be less efficient than that prediction because of their underlying 
dynamical character. Hence, it is argued that the fast electronic motion embedded 
in quantum fluctuations of atom dipoles cannot be instantaneously followed by free 
charges because of inertia. According to that heuristic argument, it is oftenly as- 
sumed that U v aw(R) is not screened at all. Notice that such argument is close to 
the one presented previously which suggests a breakdown of exponential screening 
in the quantum case. 



III. SOME EXACT RESULTS WITHIN THE SCREENED CLUSTER 
REPRESENTATION 

A. Fundamental issues and the hydrogen quantum plasma 

In the previous introductory part of that lecture, various fundamental issues have 
been raised, namely 

• Nature of screening in quantum systems 

• Introduction of recombined entities in the many-body problem 

• Modification of van der Waals interactions due to free charges 

According to the preliminary analysis and discussions, it clearly appeared that a 
proper and reliable account of the relevant mechanisms requires the introduction of 
an elementary description of matter in terms of a quantum Coulomb system made 
with point nuclei and point electrons. For the sake of simplicity, and also because 
of its wide interest for both conceptual purposes and practical applications, here we 
will restrict ourselves to the hydrogen plasma viewed as a system of quantum point 
particles which are either protons or electrons, interacting via the instantaneous 
Coulomb potential. Protons and electrons have respective charges, masses, and 
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spins, e p = e and e e = — e, m p and m e , cr p = <7 e = 1/2. In the present non- 
relativistic limit, the corresponding Hamiltonian for N = iVp + N e particles reads 



with v (r) = l/r, (a, = p,e) the species of the ith particle, and A; the Laplacian 
with respect to its position Xj. In the following, we present several exact results 
for the equilibrium properties of that system, in relation with above fundamental 
issues. 

First, in Section IIII Bl we recall various rigorous results about the existence of 
the thermodynamic limit, as well as about the partially ionized atomic regime. 
Then, in Section IIII CI we introduce the path integral representation, in the simple 
case of a single particle submitted to an external potential for the sake of pedagogy. 
Application of that tool to the many-body quantum system leads to the introduction 
of an equivalent classical system made with extended objects called loops, which 
interact via a two-body potential equal to some average along their shapes of the 
genuine two-body particle interaction. The equilibrium quantities of the quantum 
plasma can then be represented by Mayer-like diagrammatical series in the loops 
world. In Section IIII Dl we describe the resummation machinery which exactly 
transforms the sum of divergent Mayer graphs into a sum of convergent graphs with 
a similar structure, the so-called screened cluster representation (SCR). Long-range 
Coulomb divergences are removed via systematic chain resummations, which amount 
to introduce a screened effective potential between loops. Further reorganizations 
lead to the introduction of particle clusters with finite statistical weights, which 
include the contributions of familiar chemical species. Exact asymptotic expansions 
in the partially ionized atomic regime can be derived within that SCR, for both the 
equation of state and particle correlations, as decribed in Section IIII El Eventually, 
we summarize m Section IfflFl the main answers to the fundamental issues quoted 
above, which can be inferred from the present analysis of the hydrogen plasma. 
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B. Rigorous results 



1. Thermodynamic limit 

Let us consider the hydrogen plasma in the framework of the grand-canonical 
ensemble. The system is enclosed in a box with volume A, in contact with a ther- 
mostat at temperature T and a reservoir of particles that fixes the chemical poten- 
tials equal to \i v and p e for protons and electrons respectively. The corresponding 
grand-partition function reads 

H A = Tr A exp(-/3(F iVpjJVe -/i p iV p -/i e iV e )) . (III.2) 

The trace Tr A is taken over a complete basis of N-body wavef unctions, which are 
symmetrized according to Fermi statistics and satisfy Dirichlet boundary conditions 



at the surface of the box, while particle numbers iVp and N e vary from to oo 79 |. 



The thermodynamic limit (TL) is defined by the limit A — > oo at fixed tempera- 
ture and fixed chemical potentials. In that limit, Lieb and Lebowitz [30] proved the 
existence of a well-defined bulk equilibrium state with the right extensive properties 
expected from macroscopic thermodynamics. That state is overall neutral, namely 
the average particle densities p p> A = (N p )/A and p e) A = (N e )/A become identical in 
the TL, 

limp PjA = limp ei A = p . (III. 3) 

TL TL ' 

That common particle density depends only on temperature T and on the mean 
chemical potential 

„ = , (III.4) 

while the difference v = (// e — /i p )/2 is not relevant. Pressure P defined by the 
standard formula 

P = fc s Tlim^^ (III.5) 

TL A 

is also a function of the two parameters T and //, as well as all the other intensive 
thermodynamical quantities. Moreover, bulk properties no longer depend on the 
choosen boundary conditions, and all statistical ensembles become equivalent in 
the TL. In particular, grand-canonical expression ( 1III.5j) for pressure P(T, //) does 
coincide with its canonical counterpart P(T, p) evaluated at density p = p(T,p). 
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Lieb and Lebowitz proof involves sophisticated mathematical analysis which will 
not be detailed here, of course! However, it is instructive to present the main 
ideas underlying the derivations. In fact, the central difficulties are related to the 
singular behaviours of the Coulomb potential at both short and large distances, 
as already quoted in Section III A[ The short-distance collapse is avoided thanks 
to a combination of the Heisenberg uncertainty principle and of Fermi statistics. 
According to that principle, the 1/r- singularity is smeared out, so the groundstate 
energy of A" particles is finite, instead of — oo when opposite charges are present. 
Fermi statistics prevent the modulus of negative groundstate energies to grow faster 
than N. At a mathematical level, this is rephrased in the celebrated H-stability 
theorem 31], 

H Np , Nc > -B (jV p + N e ) (III.6) 

where B is a strictly positive constant. The large-distance explosion is prevented 
by screening mechanisms, the mathematical formulation of which is embedded in 
the cheese theorem. That theorem relies on a suitable partition of space in neutral 
spheres surrounding the charges, which applies for most probable configurations. 
Indeed, such configurations are overall neutral in the bulk while excess charges 
are repelled on the boundaries, as suggested by macroscopic electrostatics. The 
harmonicity of Coulomb potential then plays a crucial role, since according to Gauss 
theorem, two non-overlapping neutral spheres do not interact. This implies that 
dangerous long-range contributions do not intervene in the TL. Both cheese and 
H-stability theorems are the key ingredients of the proof. 



2. Partially ionized atomic regime 

In Section III Dl we showed how, according to phenomenological Saha theory, 
a finite fraction of protons and electrons might recombine into hydrogen atoms, 
forming a partially ionized atomic gas. Also, we argued that such a regime should be 
attained when both temperature and density are sufficiently low so that individual 
atoms in their groundstate can form, namely kT <C \Er\ and a ^> a# with a = 
(3/(47rp)) 1 ' 3 the mean inter-particle distance. 

The above Saha prediction has been rigorously proved within a suitable scaling 
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limit where both temperature and density vanish in a related way. That scaling 
limit is introduced in the framework of the grand-canonical ensemble : temperature 
T goes to zero while average chemical potential p approaches the value Ejj with a 
definite slope, according to parametrization 32I 

li = E H + fc s T[ln( 7 ) + ln((m/M) 3/4 /4)] (III.7) 

with fixed dimensionless parameter 7 and M = m p + m e . Then density p indeed 
vanishes with T, and it behaves as 

p = p* 7 (l + 2) (1 + 0{e~*)) , (III.8) 

with c a postive constant and temperature-dependent density p* given by formula 
( 1H.25|) . The corresponding asymptotic behaviour of pressure P reads 



PP = p* 7 (2 + \) (1 + 0{e-^)) . (III.9) 

Thus, if we discard the exponentially small corrections embedded in the 0(e~ c/3 )- 
terms, the corresponding leading behaviours of density and pressure can be reinter- 
preted in terms of ideal densities of protons, electrons and atoms, identified to 

pf ] = P ( : d) = p*7 an. 10) 

and 

pl?=P*y- (111.11) 



Such expressions do coincide with their Saha counterparts ( 111.231) and (1II.24|) . as it 



can easily checked by eliminating 7 in favor of p. Also, the pressure indeed reduces 
to that of an ideal mixture of protons, electrons and atoms, and it eventually reads 

/3Psaha = P* (p/p* + y/l + 2p/p* - l) . (III. 12) 

That rigorous derivation constitute another tour de force in mathematical physics, 
which has been implemented through successive works 32j, |33] after the pionneering 
paper about the purely atomic limit by Fefferman {34]. Below, we sketch various 
estimations and arguments which provide a simple physical understanding of the 
result. 
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First, the above scaling limit defines the proper energy-entropy balance which 
ensures the emergence of the contributions of ionized protons, ionized electrons and 
hydrogen atoms in grand-partition function (IIII.2j) . That balance results from the 
competition between the fugacity factors exp(/3/i) which decay exponentially fast as 
exp((3E}i), and the groundstate Boltzmann factors which may explode exponentially 
fast. The relevant ideal contributions of ionized protons and ionized electrons are 
merely given by terms (N p = l,N e = 0) and (N p = 0,N e = 1), and they are of 
order exp(/3i?H)- The corresponding atomic contribution, easily picked out in term 
(iVp = 1, N c = 1), is also of order exp(/3_E H ) since it involves the product of entropy 
factor exp(2/3/x) times the atomic Boltzmann factor exp(— PEn). Thanks to remark- 
able inequalities [80 1 satisfied by the groundstate energies of Coulomb Hamiltonians 
Hn p ,n c for iVp + N c > 3, all the other contributions of more complex recombined 
entities decay faster than exp((3En). For instance, the molecular contribution em- 
bedded in term (N p = 2, N e = 2) is of order exp(4/3£ , H) exp(— f3Eu 2 ), and it decays 



faster than exp(/3_E H ) by virtue of inequality [81] 



3E U < E H2 . (III. 13) 

Once, the contributions of molecules and other complex recombined entities have 
been discarded, it remains to estimate exchange and interaction contributions of 
ionized charges and atoms. In fact, since density vanishes exponentially fast, typ- 
ically as exp((3Eu), the mean-interparticle distance a = (3/(47rp)) 1//3 grows much 
faster, when T — > 0, than all thermal de Broglie wavelengths A p , A e , A at of ionized 
charges and atoms. Also, the average interactions between those entities, which are 
proportional to inverse powers of a, vanish exponentially fast when T — > 0, so they 
become negligible compared to thermal energy fcgT. Thus exchange and interaction 
contributions can be also dropped out in the scaling limit of interest. According to 
the previous simple estimations and rough arguments, it is not surprising that the 
system behaves as an ideal Maxwell-Boltzmann mixture of ionized protons, ionized 
electrons and atoms, in the considered scaling limit. 

Let us define the so-called Saha regime by low but finite temperatures T, and 
very low but finite densities p of order exp(/3i?H)- According to the previous rigor- 
ous result, hydrogen should behave as a partially ionized atomic gas, weakly coupled 

27 



and weakly degenerate, including also small fractions of other complex entities, like 
molecules H2 or ions and H~ for instance. As illustrated further, it is particu- 
larly interesting, at both practical and conceptual levels, to derive the corresponding 
corrections to the ideal terms obtained within Saha theory. For that purpose, it is 
convenient to combine the path integral representation to diagrammatical tools. 
Notice that the mathematical techniques involved in previous proof do not pro- 
vide explicit expressions for the exponentially small corrections beyond ideal Saha 
formulas. 



C. Path integral representation 

1. Feynman-Kac formula 

For the sake of pedagogy, we introduce the path integral representation in the 
simple case of a single particle with mass m submitted to a potential V(r). Its 
Hamiltonian reads 

h 2 

H = -— A + V(r). (111.14) 
2m 

The corresponding density matrix at a given temperature T, namely the matrix 



element of Gibbs operator exp(—/3H), is exactly given by Feynman-Kac formula 35 



x exp[-/3 / d S y((l-s) r a + sr b + A $(s))} , (111.15) 
Jo 

with thermal de Broglie wavelength A = ((3h 2 /m) 1//2 . In the r.h.s. of flIH.15j) . £(s) is 
a dimensionless Brownian bridge which starts from the origin at dimensionless time 
s = and comes back at the origin at dimensionless time s = 1, i.e. £(0) = £(1) = 
0. Functional measure is the normalized Gaussian Wiener measure which 

characterizes the Brownian process, and it is entirely defined by its covariance 

J V{£) £ M (s) Ut) = V inf M) (1 - sup( S ,t)) • (IIL16) 

The corresponding functional integration is performed over all Brownian bridges 
£(s). Representation flIII.151) is the proper mathematical formulation of genuine 
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Feynman's idea, which amounts to express the density matrix as a sum over all 
possible paths going from r a to rb in a time (3h, of weighting factors exp(—S/h) 
where S is the classical action of a given path computed in potential —V. Here, 
such paths are parametrized according to 

u ab (sf3K) = (1 - s) r a + s r b + A £(s) , (111.17) 

where (1 — s) r a + s r b describes the straight uniform path connecting r a to r b (see 
FigHl). Also, weighting factor exp(— S/K) is splitted into the product of three terms. 
The first term, which arises from the kinetic energy of the straight uniform path, 
reduces to the Gaussian prefactor in front of the functional integral. The second 
term, associated with the kinetic contribution of the Brownian part of the path, is a 
Gaussian functional of embedded in Wiener measure The third and last 

term is rewritten as the Boltzmann-like factor associated with time average 

ds V (w ab (s/3/i)) (III.18) 







of potential V along the considered path u ab . We stress that, independently of 
the rather poetic introduction of path integrals by Feynman 39(, representation 
(IIII.15I) can be derived in a straightforward way by starting from the obvious identity 
exp(— f3H) = [exp(— f3H/N)] N combined with a suitable insertion of (iV — 1) closure 



relations in position-space (see e.g. Ref. jf}]). Feynman-Kac formula (jlll.15]) then 
follows by taking the limit N — > oo, as it has been proved for a wide class of 
potentials j^]. 

Feynman-Kac (FK) representation ( 1III.15I) perfectly illustrates the intrusion of 
dynamical features in equilibrium static quantities for quantum systems, which 
has been quoted in Section III CI Let us consider the diagonal density matrix 
(r a | exp(—(3H)\r a ). Because of the non-commutativity of the kinetic and potential 
parts of H, that matrix element does not reduce to its classical counterpart 

(2ttA 2 )!>/ 2 ' 1 ' 

so it is not entirely determined by V(r a ). In fact, according to formula f lIII.15j) spec- 
ified to r b = r a , (r a | exp(— f3H)\r a ) now depends on the potential landscape in some 
neigbourhood of r a with size A, which is explored in time-average (1III.18I) thanks to 
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FIG. 4: A brownian path t*; a b(s/3?i). The dashed straight line is the uniform path connect- 
ing r a to r b . 

Brownian motion. Therefore, (r a | exp(— f3H)\r a ) appears to be indeed generated by 
a dynamical process. Not surprisingly, particle mass m controls the importance of 
the corresponding dynamical effects. In particular, in the limit of an infinitely heavy 
particle m — > oo, A = {(3h 2 /m) 1 / 2 vanishes and (r a | exp(— (3H)\r a ) obviously tends 
to its classical counterpart (1III.I9I) : dynamical effects do not intervene anymore in 
the potential contribution which takes its purely static form j^. Notice that the 
present dynamical considerations are the manifestation, in the framework of path 
integrals, of the Heisenberg uncertainty principle which prevents the particle to stay 
at r a . Hence, Brownian paths can be interpreted as describing intrinsic quantum 
fluctuations of position. 

Remarkably, FK representation ( 1III.1 51) involves only classical objects and c- 
numbers, so the operatorial structure of quantum mechanics is, in some sense, erased. 
That feature turns out to be particularly useful in the framework of the many-body 
problem, as described further. However, the intrinsic complexity of quantum me- 
chanics is now hidden in the functional integration over all Brownian bridges, which 
remains a formidable task. In fact, explicit calculations can be performed in a 
few number of reviewed in Ref. jssj]. Also, aymptotic Wigner-Kirkwood 

expansions of (r a | exp(— f3H) |r a ) around the classical formula (1IH.I9|) . can be de- 
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rived for situations where A becomes small compared to the characteristic length 
[|W(r a )|/|V(r a )|] _1 of variation of V(r) in the neighbourhood of r a . Such situations 
occur for very heavy particles or for high temperatures, and also at large distances 
for potentials which decay as power laws at infinity In the opposite limit where A 
diverges, direct estimations of FK functional integral become rather cumbersome. 
Notice that, if H has a single isolated boundstate with energy Eq and wavefunction 
ipo, in the zero-temperature limit, the asymptotic behaviour of (r a | exp(— /3if)|r a ) is 
merely extracted from its spectral representation, i.e. 

(r a |exp(-/3#)|r a ) ~ |^ (r a )| 2 exp(-/3£ ) when T^O. (111.20) 



As argued in Ref. 43], the relevant paths which provide the low-temperature be- 



haviour (IIII.20j) occupy a small piece of the whole functional phase space, and they 



are quite different from the typical paths with divergent size A. Consequently, an 
exact direct estimation of their contribution remains an open problem in general. 

2. Gas of loops 

Let us come back now to the hydrogen plasma in the framework of the grand- 
canonical ensemble described in Section IIIIBI The trace in the grand-partition 
function ( 1III.2j) can be expressed in the basis of positions and spins, where a given 
state is the antisymmetrized Slater product of one-body states |x a^)- This provides 
a sum of diagonal and off-diagonal matrix elements of exp(— /3i?jv p ,AT e )- An example 
of such matrix element with (iVp = 3, iV e = 4) is 

(R 1 R3R 2 r 2 r3r 1 r4|exp(-/3if3,4)|RiR2R3rir 2 r 3 r 4 ) , (111.21) 

where the positions of two protons are exchanged, as well as those of three electrons. 
Contributions of spins are factored out in simple degeneracy factors because the 
Coulomb Hamiltonian ^ e does not depend on the spins. For matrix element 
( 1IH.21|) . that multipliyng degeneracy factor is 2 4 , because the spin-states of the 



exchanged particles are necessarily identical. 

The FK representation for each of the matrix elements of exp(—/3HN p ,N c ) takes 
a form similar to formula (IIII.15I) . with N p protonic paths and N e electronic 
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FIG. 5: A set of four loops constructed from matrix element (|III.21j) . 



paths uj( e \ as well as a Boltzmann-like factor associated with the time average of the 
potential part Vn p ,n c of Hn p< n c . The paths associated with matrix element (1IH.21|) 
are drawn in Fig. We see that all those paths can be collected in loops. In 
fact, that property holds for any matrix element of exp(— f3H^ p jv c ), because any 
permutation can always be decomposed as a product of cyclic permutations. A loop 
£ is constructed by collecting q paths associated with q particles exchanged in a 
cyclic permutation. Accordingly, C is characterized by its position X, which can be 
arbitrarily chosen among the extremities of paths lo, and several internal degrees of 
freedom which are particle species (a = p,e), number q of exchanged particles, and 
shape X a if obtained as the union of the q paths uj. It turns out that r](s) is itself 
a Brownian bridge with flight time q, i.e. f](0) = 77(g) = 0, distributed with the 
corresponding Wiener measure T>(rj). 

In the FK representation, the time-average of potential Vjv p ,N e can De obviously 
rewritten as a sum of two-body interactions V between loops, plus a sum of loop 
self-energies 14, namely 



where loops associated with the considered matrix element are labelled as with 
index i running from 1 to their total number N. For instance, four loops can be 
identified in the FK representation of matrix element ( 1III.21I) . Two-body potential 
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V(£j, £j) between loops and Cj reduces to a time-average along their respective 
shapes of the genuine two-body particle interaction v(\Xi + \ ai r]i(s) — ~Kj — \ aj r]j(t)\) 
evaluated at times which differ by an integer value. Self-energy U{Cj) for loop £j is 
given by a similar average along its own shape of v(\\ ai rji(s) — X ai r}i(t)\) evaluated 
at times which differ by a non-zero integer value, with a prefactor 1/2 which avoids 
double counting of genuine interactions between two exchanged particles. 

At this stage, grand-partition function (1III.2j) is rewritten as a sum of Boltzmann- 
like factors associated with energies ( 1III.22I) multiplied by combinatorial factors and 
particle fugacities z a j (27rA 2 ) 3 / 2 with z a = exp(/3/z Q ), which have to be integrated 
over positions and shapes of the involved loops. It turns out that the whole sum can 
be rewritten as the grand-partition function of a classical gas of undistinguishable 
loops with suitable activities z(C), namely 



Phase space measure dC in the world of loops, involves discrete summations over 
species index a and exchanged-particle number q, spatial integration over position 
X inside A, and functional integration over shapes rj with Wiener measure T>(rj) 
restricted to shapes such that X + XaV( s ) remains inside A. Loop fugacity reads 



the structure of which is easily intepreted as follows. Factor (— is the signature 
of a cyclic permutation of q objects. In front of the obvious particle-activity con- 
tribution z^, factor 2 is the number of configurations of exchanged-particles spins 
which are all identical. Factor q in front of (27rgA 2 ) 3 / 2 is related to the q possible 
choices of position X among that of the q exchanged particles. The other factor 
q in front of A 2 , arises from the absorption of the remaining (q — 1) integrations 
over particle positions together with the q functional integrations over Brownian 
bridges into the single measure T>(rj). Notice that q\ 2 a is nothing but the square of 
de Broglie wavelength for inverse temperature q/3. 

We stress that identity ( 1III.23I) proceeds from remarkable combinatorial prop- 
erties. For instance, in S^ oop \ there are various contributions in term N = A, 




(111.23) 



z{C) = (-1) 



g(2vrgA 2 ) 3 / 2 



exp(-pU(C)) , 



(111.24) 
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which are identical to that of matrix element ( lIII.2ip in the Slater expansion of Ha- 
Such contributions arise from the explicitation of the 4 phase-space measures dCi 
(i = 1,2,3,4) where 2 protonic loops and 2 electronic loops, carrying respectively 2 
or 1 protons and 3 or 1 electrons, have to be chosen among 4 labelled loops. Also, 
there are several matrix elements in the Slater expansion of 2a, involving 3 protons 
and 4 electrons, which provide identical contributions to that of matrix element 
(IIII.21I) . Thanks to a rather fortunate arrangment between both counting factors, 



the full respective contributions in both H A loop ^ and 2a are indeed identical! Accord- 
ing to those considerations, and also because of its synthetical form, identity (1IH.23|) 
is oftenly called the magic formula. 



Historically, Ginibre [44| was the first to introduce the notion of loops when 
studying the convergence of Mayer series for quantum gases with short-range forces. 
However, he did not write explicitely formula ( 1III.23I) . which has been derived later 



by Cornu [45(|. Recently, Martin [46| proposed an elegant and shorter derivation 
of that formula. Notice that identity (1III.23I) is valid for any kind of two-body 
interactions, and an arbitrary number of species with Fermi or Bose statistics. For a 
bosonic species, factor ( — is merely replaced by 1 in loop fugacity flIII.24j) . while 
all the other factors are unchanged. Magic formula (1111.230 has been applied to a 



Bose gas with short-range interactions for studying Bose-Einstein condensation [47 ]. 



Of course, and as quoted in the simple case of a single particle in an external 
potential, the intrinsic difficulty of quantum mechanics is now hidden in the func- 
tional integrations over loop shapes, so an exact calculation of loop grand-partition 
function remains far beyond human abilities... Nevertheless, magic formula ( 1III.23j) 
is quite useful because standard tools of classical statistical mechanics can be ap- 
plied, as well as various transformations relying on simple properties of classical 
Boltzmann factors. 



D. Screened cluster representation 

From now on, we assume that the thermodynamic limit has been taken once 
for all. Since, according to Lieb and Lebowitz theorems, boundary effects do not 
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intervene anymore in the TL, we proceed to formal calculations in the infinite system 
where some quantities are infinite. Nonetheless, after suitable resummations and 
reorganizations, such divergences are removed, and it is quite reasonable to believe 
that the final results are indeed relevant, as far as they are finite, and they do 
describe the considered bulk quantities. That strategy is adopted for the calculation 
of particle densities and particle correlations. Also we work in the world of loops. 
According to magic formula (IIH.23I) . those particle quantities are merely related to 



their loop counterparts. For instance, proton density p p is given in terms of loop 
density p(£ a ) by 

CO „ 

Pp = E / <?a P(£a) , (IIL25) 

where loop species index is a a = p. 
1. Mayer graphs 

The structure of loop grand-partition 23[) is identical to that of an ordinary 
classical system made of point particles with two-body interactions. Therefore, 
equilibrium quantities of loops can be represented by Mayer-like diagrammatical 
series, where points are replaced by loops. This provides Mayer series for particle 
quantities, like proton density p p which reads 

Pp = E ^7n) E / *W / II d£ * [II /] a ■ ( IIL26 ) 

Q ^ ' q 3 ,=X J J i=l 



Each graph Q is constructed according to the standard Mayer rules [48|, |49|. It 
is made with (n + 1) loops £j, i = 0, ...,n and £ = £&■ Two loops i and j are 
connected at most by a single bond 

^=exp(-/3V(A,£,))-l, (111.27) 

and [Yl f]g denotes the product of such bonds. Also, graph Q is simply connected, 
namely it cannot be separated into two parts which are not connected by at least 
one bond /. Symmetry factor S(Q) is the number of permutations of black loops Li 
with i > 1, which leave Ylg f unchanged. Each loop is weighted by its fugacity z(Ci). 
The contribution of graph Q is obtained by integrating over all degrees of freedom 
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Co 



FIG. 6: A graph Q made with root loop C a and two field loops C\ and £2- 

of black loops embedded in measure d£«, while the position and species index of 
root loop £ a are fixed and only its shape and particle number are integrated over. 
Eventually, Ylg 1S performed over all topologically different unlabelled graphs Q. An 
example of graph Q is shown in Fig. El while its contribution reads 

lit [ 9a z{C & ) I d£id£ 2 zi&WCJUfa ■ (IIL28) 

The loop-loop interaction is long ranged, like the genuine Coulomb potential 
itself, as illustrated by the large distance behaviour 

V(A, Cj) ~ qj L aSj( ^ ] when |X, - Xj| 00 . (111.29) 

Indeed, at large distances, loops can be shrinked to point charges, and further correc- 
tions to the monopolar form (lIII.29j) can be expanded in multipolar power series of 



A Qi /|Xj — Xj| and A a ./|Xj — Xj|. Like in the case of purely Coulombic interactions, 
the long-range nature of V(Ct,Cj) induces divergences as explained in Section Hi Al 
for purely Coulombic interactions. Here, in Mayer series 26j) . such divergences 
pollute every graph because Mayer bonds behave as — (3V(£i, Cj) at large distances. 

In addition to previous long-range pollution, another drawback of Mayer series 
relies on the difficulty to identify immediately the contributions of a given chemical 
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Resummation Machinery 

Expand all Mayer bonds / in Taylor series with respect to (— f3V) 

Resum all convolution chains of bonds (— /3V) connecting Li to Cj 
in terms of an effective potential e ai e aj cp(Ci, Cj) 

Introduce Boltzmann factors exp(—{3e ai e aj <j)(£i,£j)) and renormalized 
fugacities z^(Ci) via suitable reexponentiations 

Identify sets of loops 6c = {C-i} such that all the effective interactions between loops 
inside 6c are present in the product of Boltzmann factors ex.p(—(3e ai e a:j (fr(Ci,Cj)) 

Identify particle clusters C(N p , N e ) in a given set 6c 
by expliciting dCi as sums over species index and particle numbers 

Introduce statistical weights ZT(C) for particle clusters and bonds 
via suitable truncations of Boltzmann factors exp(— f3e a -e a (j)(Ci, Cj)) 

• * 



FIG. 7: Main steps of the resummatiom machinery. 

species. For instance, if we are interested in contribution of ion H^" to series (IIII.26I) . 
we can pick out in contribution (1111.281) of graph drawn in Fig. [BJ the term (q a = 
l,a a = p ; qi = l,«i = p ; q<i = 1,«2 = e), which involves 2 protons and 1 
electron. However, because of the absence of bond /12 in that graph, a proton- 
electron interaction is missing. The full interaction between the two protons and 
the single electron is recovered by summing contributions from all graphs Q made 
with three loops. 

2. Resummation machinery 

The resummation machinery amounts to perform an exact transformation of the 
whole series (IIII.26j) . which takes care of drawbacks of Mayer grahs quoted above. 
The main steps are sketched in Fig. [71 Of course, all the transformations involve 
tricky counting calculations based on combinatorial identities, which are detailed in 
Ref. [50|. The final result can be expressed as the Screened Cluster Representation 
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c«(i,i) 



Ci(l,l) 



FIG. 8: A graph G made with root cluster C a (l, 1) and one field cluster Ci(l, 1). Each 
cluster contains one proton and one electron denoted by symbols p and e respectively. 

(SCR), 

pv=Y<s{G)J dc * & ^(cy j n ^(co [n ^} q > ( nL3 °) 

where graphs G have the same topological structure as usual Mayer graphs. Now, 
ordinary points are replaced by particle clusters Cj, weighted by factors ZT(d). Two 
clusters are connected by at most one bond which can be either — /3 2 $ 2 /2! 
or — /3 3 $ 3 /3! with $ the total effective potential between both clusters. Graphs G 
are simply connected, and symmetry factor S(G) is computed as usual. Phase-space 
measure dCj reduces to integrations over positions and shapes of loops belonging 
to set £c associated with cluster Cj, whereas in dC a position of root proton is not 
integrated over. Eventually, YIg * s performed over all topologically different graphs 
with some simple restrictions, in part avoiding double-counting of genuine particle 



interactions 50[- In Fig. EJ we give an example of graph G, made with two clusters 



C a = C a (l, 1) and C\ = Ci(l, 1). Its contribution reads 



dC7 a Zj(C a ) / dd Z T JC X ) ^^C!) ^ (nL31) 



since S(G) = 1 and q a = 1. Notice that part of graph Q shown in Fig. [6] contributes 
to previous graph G under the action of the resummation machinery. 

Remarkably, not only the genuine structure of Mayer graphs is conserved through 
the resummation machinery, but both statistical weights Zj(C) and bonds de- 
pend on the sole effective potential 0. Thanks to the sufficient fast decay of (p, 
all large-distances divergences are removed in every graph G, which does provide 
a finite contribution. The presence of complete interactions inside a given cluster 
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C(N p , N e ) ensures that chemical species made with N p protons and N e electrons do 
emerge at low densities. Both features are detailed below. Notice that the physical 
ideas underlying the SCR, are close to those involved in the construction of the 



so-called ACTEX approach by Rogers 5JJ. Also, contrarily to Feynman graphs in 



volved in standard many-body perturbation theory [20|, graphs G in SCR account 
for non-perturbative effects in the Coulomb potential which are essential for a proper 
description of recombination. 

3. Screened effective potential 

Effective potential e aa eo, b 0(£ a , £ b ) is defined by the sum of all convolution chains 
of bonds (—/5V) connecting £ a to £ b , namely 

oo „ n 

e a ^e ab <p(C a , £ b ) = V(£ a , £ b ) + ^(-1)"/?" / \{ d£, z{L { ) V(£ i _ 1) A) V(£ n , £ b ) , 

n=l J i=l 

(111.32) 

with £ = £ a . That potential accounts for collective effects, because it involves 
contributions of an infinite number of loops. Interestingly, chain sum flIII.32j) can 



be reinterpreted as the complete perturbative expansion of the solution of integral 
equation 

e Qa e ab 0(£ a , £ b ) = V(£ a , £ b ) - /? J d£ z{C) e Qa e a 0(£ a , £)V(£, £ b ) . (111.33) 

Therefore, eo, a e ab 0(£ a , £ b ) is the total potential acting on loop £ b , created by loop 
£ a plus a cloud of loops £ distributed with density — /3z(£) e ai ,e a (j)(C a , £). Its 
physical content is then quite similar to that of the Debye potential introduced in 
Section III Bl In the classical limit fi — > 0, loops shrink to point charges so equation 
(1III.33I) is nothing but the integral version of linearized Poisson-Boltzmann equation 
( 111.91) with e w = 1. Thus, the effective potential e Qa e ab 0(£ a , £ b ) can be viewed as 
the quantum analogue of the Debye potential [521 ] . 



Similarly to what happens in the classical case, collective effects embedded in the 
effective potential are expected to screen the bare loop- loop interaction. In fact, the 
static part g a e Qa g b eQ, b /|X a — X b | of V(£ a , £ b ), which does not depend on loop shapes, 
is exponentially screened fs3 1 - On the contrary, the multipolar parts of V(£ a ,£ b ) 
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which depend on loop shapes, are not exponentially screened. The resulting effective 



potential decays only as a dipolar interaction at large distances, namely 52] 



0(£ a , 4>) ~ 4^%T when |X a -X b |^oo, (111.34) 
|X a - X h \ 6 

where -4a a a b (f7a, 1%) is an amplitude which depends only on loop species and shapes. 
Consequently, bonds = — /3<& are at the border line for integrability at large 
distances. However, if loop shapes are first integrated over, all dangerous l/R 3 - 
terms vanish. Thus, within that prescription for the order of integrations involved 
in measures dCj, each graph G in SCR flIII.30j) indeed provides a finite contribution. 



The slow algebraic decay of can be understood as follows. As mentionned in Sec- 
tion III I C\ paths defining loop shapes account for intrinsic quantum fluctuations of 
particle positions. Such fluctuations, which are dynamical in origin, generate 1 / R 3 - 
dipolar interactions. In agreement with heuristic findings argued in Section III CI 
such dynamical interactions cannot be perfectly screened trough collective effects, 
so it remains a 1/ i? 3 -tail in the large- distance behaviour of 0. 

4- Emergence of chemical species 

As quoted in Section lHD} familiar chemical species may emerge, strictly speaking, 
in a zero-density limit only. Not surprisingly, in such a limit, effective potential 
e aa e a) 6iC a ,, £b) reduces to the bare loop-loop interaction V(£ a ,£b) as shown in 
Ref. 52j. Let us consider the contribution of the simple graph G made with the 
single root cluster C a (iVp, N e ). Functional integrations over loop shapes of products 
of Boltzmann factors exp(— {3V(£i, Cj)) embedded in Zj(C a ), can be expressed in 
terms of matrix elements of Gibbs operators associated with Coulomb Hamiltonians, 
by applying backwards FK formula. Then, the contribution of the considered graph 



G is proportional to 



(2w\ 2 r - W 2 



Z(N P , AQ = ^""y Tr [eM-PH Np , Ne ) - ...] (111.35) 

with \% n c = /3h 2 /(N p m p + N e m e ). Terms ... left over in the r.h.s. of definition 
( 1III.35j) indicate a suitable truncation of Gibbs operator exp(— (3HN pt N c ), which en- 



sures the finiteness of the trace. Thus, we see that SCR leads to a natural definition 
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of the partition function Z(N p , N e ) for N p protons and N e electrons in the vacuum. 
That partition function includes the contribution of possible recombined species in 
their groundstate, as well as those of thermal excitations. 

At finite densities, Zj(C a ) accounts for many-body effects on cluster C a , like the 
broadening and shift of energy levels of recombined species made with N p protons 
and N e electrons. Also, graphs G with several clusters describe interactions between 
such clusters which are screened by ionized charges. For instance graph shown in 
Fig. [HI involves contributions from interactions between two atoms. In fact, all 



phenomena at work can be identified in specific graphs, as argued in Ref. 50 ]. 



E. Asymptotic expansions in the Saha regime 

1. Equation of state 

Now, we consider again the Saha regime, where hydrogen behaves as a partially 
ionized atomic gas. We start with a grand-canonical description, where the given 
thermodynamic parameters are p and T, and common particle density p = p p = p c 
is a function p(p, T). Making the variable change (p, T) — > (7, T) defined through 
parametrization (1111.71) . we see that density can also be viewed as a function of 7 
and T. Then, we proceed to the expansion of p(7, T) / p* when T — > at fixed 7, 
by investigating the corresponding behaviour of graphs G in SCR (IIII.30j) . Roughly 



speaking, that behaviour results from the competition between three kind of contri- 
butions, which vanish or explode exponentially fast when T — > 0, namely 

(i) 7^ exp((3NEft) with N = N p + N e arising from entropy factor exp((3pN) in 

(ii) exp(— /3E^ N J arising from a recombined entity with groundstate energy E^' Nc 
in Zl{C) 

(iii) 7 P//2 exp((3pEu/2) with p relative integer, arising from interactions screened over 
Debye length A D = {4nf3e^ip p id) + p^))" 1 / 2 with p p id) = p£ d) = p* 7 in bonds 



between charged clusters 



84] 
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Previous analysis provides the Scaled Low Temperature expansion of pi^y, T)/p* 
where leading terms 



7 + 



r 



(111.36) 



are indeed those predicted by Saha theory. Each correction reduces to some power 
of 7 multiplied by a temperature-dependent function which decays exponentially 
fast when T — > 0, in agreement with the rigorous behaviour (1III.8H presented in 
Section UTTBl 

The SLT expansion of (3P/p* is readily obtained from that of p(7,T)/p*, by 
integrating identity 



d/3P 
(97 



09,7) 



2p 

7 



(111.37) 



whichs follows from the standard thermodynamical relation between p, P and p. A 
further elimination of 7 between both SLT expansions of /3P(7, T) / p* and p(7, T) / p* 
provides the SLT expansion of the equation of state, i.e. the asymptotic expansion 
of pressure in units of p*k^T when T — > at fixed ratio p/p*, 



f3P/ P * = /3P S aha/p* +J2 b k(p/P*)MP) ■ 



(111.38) 



fc=i 



The leading term is indeed given by Saha formula (IIII12I) . Coefficients bk{p/p*) 
are algebraic functions of ratio p/p*, while temperature-dependent functions cek(/3) 
decay exponentially fast when T vanishes, ctk(P) ~ exp(— /36k) except for possible 
multiplicative powers of j3. Expansion ( 1III.38I) is ordered with respect to increasing 
decay rates, < 8 \ < 62 < ■■■■ The first five corrections have been computed in 



Ref. [53J (see Ref. |54j for a pedagogical presentation of the calculations). In the 
following table, we summarize their physical content, as well as the expressions and 
values of the corresponding decay rates which are merely obtained by taking the 
products of above exponential factors (i), (ii) and (iii). 



Correction (fc) 


Physical content 


8 k (in eV) 


1 


plasma polarization around ionized charges 


\E H \/2 ~ 6.8 


2 


formation of molecules, atom-atom interactions 


\3E n - E n2 \ ~ 9.1 


3 


atomic excitations, charge-charge interactions 


3|P H |/4 ~ 10.2 


4 


formation of ions, atom-charge interactions 


\2E R -Err+ \ ~ 11.0 


5 


fluctuations of plasma polarization 


\E B \ ~ 13.6 
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Remarkably, all non-ideal corrections to Saha equation of state evocated in Sec- 
tion HUB] are properly identified and ordered in exact SLT expansion ( 1III.38I) . Con- 
tributions of usual chemical species do emerge in the present double low-temperature 
and low-density limit, in agreement with heuristic findings presented in Section Hi Dl 
Their relative importance, determined by the ordering of decay rates 5%, results from 
subtle inequalities between groundstate energies Nc and E# in the vacuum. For 
instance, formation of molecules H2 dominate that of ions or H~, while contri- 
butions of more complex entitites, like H^, or H3, decay exponentially faster 
than exp(— /5|J5h|) as detailed in Ref. 53] . Thermal excitations of atoms, molecules 
and ions are accounted for in functions ot-k{0) Via the corresponding cluster parti- 
tion functions Z(l,l), Z(2, 2), Z(2,l) and Z(l,2). However, we stress that such 
contributions of recombined entities are entangled to that of their dissociation prod- 
ucts. For instance, in Z(l, 1), contributions from excited atomic states cannot be 
separated from that of ionized protons and ionized electrons. Also, those cluster 
partition functions, built with bare Coulomb Hamiltonians, are indeed finite thanks 
to the substraction of non-integrable long-range parts, which are ultimately screened 
by ionized charges. Thus, only the full contribution of those species partition func- 
tions and of their screened long-range parts, makes physical sense. For instance, 
term k = 3 involves both contributions of Z(l,l) and of screened proton-electron 
interactions. In that context, the extraction from ots(/3) of a purely atomic contri- 
bution, for instance given by the phenomenological expression ([IL2ZD of Zj cr) or by 



the Planck-Larkin formula 55J, remains arbitrary. We stress that such arbitrariness 
does not cause any trouble here, since only the full contribution embedded in a 3 ((3) 
is relevant for thermodynamics. 



Eventually, notice that ordinary virial expansions in powers of p at fixed T 56 



59], can be easily recovered from SLT expansion ( 1III.38I) by expanding coefficients 
bk{p/p*) in powers of p/p*. This has been explicitely checked up to order p 2 , by 
noting that partition functions Z(l,l), Z(2, 0) and Z(0, 2) are merely related to 



two-body quantum virial functions first introduced by Ebeling 



56). 
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2. Particle correlations 



The SCR also provides diagrammatical series for equilibrium particle correlations, 
Pa a a b (r a , r b) = pa aQb (r a ,r b ) — p^ Pn, h , where graphs are similar to those introduced 
above for particle densities |50j. Not surprisingly, the asymptotical behaviour of 
pSab( r aj r b) when R — |r a — rb| — > oo, is determined by the large-distance be- 
haviour of bonds T§. It turns out that bonds — which might give a priori 
l/i? 3 -contributions, ultimately provide short-range terms thanks to the rotational 
invariance of statistical weights of particle clusters combined to the harmonicity of 



Coulomb potential 



601 ] . On the contrary, bonds /3 2 $ 2 /2! provide l/i? 6 -contributions 



60|, all 



which do not cancel out. Thus, according to that graph by graph analysis 
particle correlations are expected to decay as 1/R 6 when R — > oo, namely 

e^r b )- ^f P) , (111.39) 

whith temperature- and density-dependent amplitudes A a ^ ah ((3, p). 

In the Saha regime, amplitudes A a&ah (f3, p) can be determined within the method 
used for deriving the SLT expansion of the pressure ( 1III.38I) . Here, it is important to 



select first the graphs which contribute to the 1 / _R 6 -tails, and afterwards to take the 
scaled low-density and low-temperature limit of the corresponding contributions. 
That order of limits ensures that collective screening effects arising from ionized 
protons and ionized electrons are indeed taken into account. Notice that Debye 
screening length Ad associated with those almost classical ionized charges diverges 
in the SLT limit. 



Within above procedure, SLT expansions of A aaab (f3, p) can be derived 60]. At 
leading order, A pp ((3, p) reduces to a quadratic form in the ideal densities (IIILlOp of 
( 1III.11I) of ionized protons and hydrogen atoms respectively, 



A PP (/3, p) = pi? pi? C at _ at (T) + p p id) p a ?C P -at(T) + p a ?p(?C at _ p (T) + ... (111.40) 
where terms left over decay exponentially faster than (p*) 2 . That structure is easily 

(T) 

interpreted by noting that each proton in correlation p pp may be either ionized or 
recombined into an atom. Then, the insertion of that leading quadratic structure 
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of amplitudes in asymptotical behaviour (1III.39|) provides a natural definition of 
effective potentials between ionized charges and atoms. In particular, the atom- 
atom effective interaction is 



Taking into account the expression of C at _ at (T) derived in Ref. 60], and definition 



( III. 301) of van der Waals potential, that effective interaction can be rewritten as 



u£% t (R) = (l " c a t^j) U vdW (R) (IIL42) 

where c at is a numerical positive constant. Thus, atom-atom van der Waals interac- 
tions are not perfectly screened by free charges, which only reduce their amplitude, 
in qualitative agreement with heuristic arguments presented in Section III Dl Re- 
markably, reduction factor (1 — c at kBT/\Eu\) does not depend on the proportion 
of free charges, since positive constant c at can be expressed in terms of the sole 
atomic spectrum |60j. The density of free charges only intervenes in Debye screen- 
ing length Ad, while formula (1III.42|) holds for R 3> Ad- Notice that, in the window 
a B <C R <C Ad, ^at-at(^) does reduce to U v aw(R) apart from exponentially small 



corrections with T arising from atomic excited states 60 1. 



F. Conclusion 

According to the exact results described here, and also to other works in the liter- 
ature, the present state of the art for the various issues about screening, recombina- 
tion and van der Waals forces is summarized below, together with some comments 
about related open problems. Part of those results are described in Ref. 6l| where 
exact results for Coulomb systems at low density are reviewed. 

• Debye exponential screening is destroyed by quantum fluctuations 

Quantum fluctuations of positions generate instantaneous electrical dipoles which 
cannot be perfectly screened by the surrounding plasma because of their dynami- 
cal character. Fluctuations of the resulting l/i? 3 -dipolar effective interactions ulti- 
mately pollute equilibrium particle correlations with 1 / i? 6 -algebraic tails. Such fluc- 
tuations are not taken into account in usual mean-field theories, like Thomas- Fermi 
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or RPA, which erroneously predict an exponential decay 62|, [63|. Previous mecha- 
nism, intrinsic to quantum mechanics, is present in any thermodynamical state. Ex- 
plicit perturbative calculations of 1/ i? 6 -tails have been performed in regimes where 



free charges behave almost classically 



60 



tances R ^> Ad- Also, it has been shown 



64|, for which such tails appear at dis- 



65], through a non-perturbative analysis, 



that the effective potential between two quantum charges immersed in a classical 
plasma does decay as l/i? 6 . Thus, all those results strongly suggest the breakdown 
of exponential screening, although a rigorous derivation is not yet available. 

Besides a proof of algebraic screening, explicit calculations of algebraic tails in 
strongly degenerate plasmas remain to be done. In particular, such calculations 
might be quite useful in condensed matter, where effective electron-electron in- 
teractions jure oftenly modelized as Dirac delta functions in order to account for 
screening 



• Recombined entities must be defined in a double zero-density and zero-temperature 
limit 

In agreement with simple findings, contributions to thermodynamical quantities of 
familiar chemical species, emerge unambiguously when both density and temper- 

uster Representation, which can be devised for any 



ature vanish. The Screened C 
mixture of nuclei and electrons [50j, is a suitable theoretical framework for evaluating 
such contributions. As it can be naively expected, when T — > 0, the leading contri- 
bution of a given recombined entity is controlled by the Boltzmann factor associated 
with its groundstate energy in the vacuum. At finite temperatures, contributions 
from thermal excitations of that entity, and from its dissociation products and their 



screened interactions, are all mixed together. This is wel 
of the equation of state of hydrogen in the Saha regime 



illustrated by the analysis 



3- 



In approaches based on the chemical picture, the choice of a suitable internal 
partition function for recombined species is a central question, which has been the 
source of many controversies since the introduction of Planck-Larkin formula (see 
e.g. Refs. 66l468l]). As far as thermodynamical quantities are concerned, the analysis 
of the various contributions in their SCR confirms the arbitrariness of such a choice. 
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More interesting, that analysis shows that, in phenomeno logical theories, the internal 
partition function of a given chemical species, should be defined simultaneously 
with the related contributions of its elementary components. The SCR itself might 
serve as an useful guide for introducing suitable modelized ingredients in chemical 
approaches. 

• Free charges reduce the amplitude of van der Waals interactions 

In the framework of the many-body problem, it is natural to define effective po- 
tentials between recombined entities from the asymptotic large- distance behaviour 
of equilibrium correlations between their elementary components. Since all particle 
correlations decay as 1/R 6 , effective interactions between neutral or charged entities 
also decay a la van der Waals as 1/R 6 . The corresponding amplitudes are con- 
trolled by quantum fluctuations of electrical dipoles, the size of which is typically 
of order either the Bohr radius Ob for recombined charges or thermal de Broglie 
wavelengths for free charges. Thus, genuine van der Waals interactions between 
atoms or molecules are not perfectly screened by free charges, as a consequence of 
the breakdown of Debye exponential screening. 

Free charges should reduce the amplitude of van der Waals interactions, as sug- 
gested by the exact calculation for hydrogen in the Saha regime. At sufficiently low 
temperatures and low densities, the corresponding renormalization factor does not 
depend on p, and is linear in T. It should be quite instructive to compare that 
first -princip les prediction to that of phenomenological approaches like Lifchitz the- 



ory 69l472| . Also, the SCR of particle correlations should provide some insights for 
reliable estimations of previous reduction factor at higher temperatures and densi- 
ties. The reduction of the amplitude of van der Waals interactions by free charges 



might have important consequences on the phase diagrams, as quoted in Ref . [73| for 



a system of biological macromolecules immersed in an highly concentrated salt. 
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[75] See for instance a one-dimensional classical gas where the attractive part of the pair- 
potential has an infinite range and a vanishingly small amplitude [4|, or lattice sys- 
tems described by Ising Hamiltonians where any two spins interact with a constant 
and vanishingly small coupling (see e.g. Ref. [5J). In those double zero-coupling and 
infinite-range limits, the mean-field treatment becomes exact. 
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electron-phonon coupling [19. |21|. Despite its obvious drawbacks, Thomas- Fermi cal- 
culation suggests that the electron-electron effective potential is shorter ranged than 
the bare Coulomb interaction. 

The fractions of atoms and molecules can be determined within a simplified mix- 
ture model of non- interacting classical particles with an internal energy E-r 2 or -^H- 
The corresponding minimization equation then reduces to the usual mass action law 
associated with dissociation reaction H2 ^ 2H. 

That standard definition of the grand-canonical ensemble may be viewed as rather 
artificial in the sense that Dirichlet boundary conditions do not describe the physical 
situation which is considered in the Gibbs construction, where particles can freely go 
in and out the volume A immersed inside a much larger container. Nevertheless, in 
the thermodynamic limit A — >• 00, boundary conditions become irrelevant for bulk 
quantities. Therefore, for our purpose it is legitimate to use definition pil.2|) . where 
particle numbers are varied by some hidden operator. 

In fact, those inequalities have been conjectured by Fefferman 34|, but they are not 



yet proved. However, all the known groundstate energies of recombined entities with 
few protons and electrons do satisfy such inequalities. 

Of course, inequality . 1 3[) also intervenes in the determination of the chemical 
composition of hydrogen via the minimization procedure described in Section III Dl 
It indeed ensures that atoms are the proeminent chemical species at sufficiently low 
densities for a given low temperature. 

Quantum corrections to classical formula (|III.19|) can be expanded in powers of h? 



according to the well-known Wigner-Kirkwood expansion [4CH42I] . That expansion 
can be easily retrieved within FK representation (jlll. 15[) . by expanding time-average 
(|III.18j) in power series of A£ and by applying Wick theorem to the calculation of the 
resulting moments of 

This can be checked by replacing all loop interactions by their static parts in chain 
sum pil.32p . The corresponding integration over loop internal degrees of freedom 
involved in each dCi provides the square of some wavenumber Ki oop . The chain sum 
then becomes a geometrical series in — /k 2 in Fourier space, and the resulting 
static part of the effective potential involves factor exp(— Ki oop |X a — Xb[). 
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[84] Since ionized protons and ionized electrons behave almost classically in the Saha 
regime, contributions of screened interactions are controlled by Debye length associ- 
ated with the corresponding density of those charges. Notice that screened potential 
e aa e ab 0(£ a , A>) precisely decays as qa,e aa q h e ah exp(-K D |X a - X b |)/|X a - X b | at dis- 
tances |X a — X b | of order Ad- The algebraic l/|X a — X b | 3 -decay of takes place at 
distances |X a — X b | much larger than Ad [531 ] . 

[85] Such a modelization of screening of Coulomb interactions between electrons, may 
be sufficient for interpreting rather complicated phenomena mainly driven by other 
mechanisms. Nevertheless, the existence of long-ranged normal attractive interactions 
between electrons might induce some observable consequences. For instance, in nor- 
mal superconductivity, interactions generated via the electron-phonon coupling are 
believed to be responsible for the formation of Cooper pairs. 
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